Many-particle interference beyond many-boson and many-fermion statistics 
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Identical particles exhibit correlations even in the absence of inter-particle interaction, due to 
the exchange (anti)symmetry of the many-particle wavefunction. Two fermions obey the Pauli 
principle and anti-bunch, whereas two bosons favor bunched, doubly occupied states. Here, we 
show that the collective interference of three or more particles leads to a much more diverse behavior 
than expected from the boson-fermion dichotomy known from quantum statistical mechanics. The 
emerging complexity of many-particle interference is tamed by a simple law for the strict suppression 
of events in the Bell multiport beam splitter. The law shows that counting events are governed by 
widely species-independent interference, such that bosons and fermions can even exhibit identical 
interference signatures, while their statistical character remains subordinate. Recent progress in the 
preparation of tailored many-particle states of bosonic and fermionic atoms promises experimental 
verification and applications in novel many-particle interferometers. 



I. INTRODUCTION 

The symmetrization postulate enforces the 
(anti)symmetrization of the bosonic (fermionic) many- 
particle wavefunction [T] and thereby severely restricts 
the set of accessible states for indistinguishable particles. 
When one postulates that each microscopic state is 
populated with equal probability [5], the resulting sta- 
tistical physics of bosons, fermions and distinguishable 
particles heavily differs, as directly observed in two-point 
correlation functions. The latter reveal bunching of 
bosons [31 3] and the opposed anti-bunching of fermions 
[5H7], which can also be directly compared in a single 
setup [S]. The differences between the species is often 
ascribed a rather universal character and said to be 
rooted in many-particle interference [7HS]- However, the 
states that are prepared in such many-body experiments 
are thermal, such that only the statistical behavior 
of bosons and fermions is probed, and no coherent 
many-particle interference crystallizes out, as we will 
explain further down. 

As a prominent example for such many-particle inter- 
ference, and in only ostensible agreement with the statis- 
tical behavior of many bosons, two single photons exhibit 
the Hong-Ou-Mandel [10] effect: Two indistinguishable 
photons that fall simultaneously onto the input modes 
of a beam splitter with reflectivity 1/2 always bunch and 
leave the setup together. Fermions behave in the opposite 
way, the Pauli principle enforces them to anti-bunch and 
to choose distinct output modes jUj . Also for more than 
two photons, bosonic effects boost the probability to find 
all particles in one output mode (T^HH]- Many-particle 
interferences thus seem to boil down to a behavior that 
is familiar from statistical physics. 

Here we falsify this popular view and show that scat- 



tering events with many particles and many modes give 
rise to much richer many-particle interference phenom- 
ena than intuitively expected from few-particle interfer- 
ence and quantum statistical mechanics. The behavior of 
bunched bosons and anti-bunched fermions, which was 
also presumed to dominate many-particle interference 
jl6| . is widely insufficient for the understanding of the 
coherent behavior of many particles, since interference 
obtrudes the overall picture. In particular, bosonic and 
fermionic interference effects are not necessarily opposed 
to each other, but fermions can experience fully destruc- 
tive interference in the very same way bosons do. Only 
when the purity of the initial state is destroyed, many- 
particle interference patterns vanish, and a smooth, fa- 
miliar, bosonic or fermionic behavior is recovered. 

We establish our results by theoretically comparing 
the scattering of distinguishable particles to bosons and 
fermions in a setup with many modes. Particles are ini- 
tially prepared in n input modes, they then scatter off a 
common potential, such that each particle ends in a co- 
herent superposition of the n output modes. The proba- 
bility for a counting event, i.e. for an event with a certain 
number of particles in each output mode, then reflects 
many-particle interference of bosons and fermions, when 
contrasted to distinguishable particles, for which usual 
combinatorial laws apply. Our system of n input- and n 
output modes can be realized with multiport beam split- 
ters for photons |17) or by an appropriate sequence of 
tunneling couplings in experiments with ultracold atoms 
in optical lattices ^TSj [19] . The distinguishability of the 
particles can be achieved by misaligning the path lengths 
for photons, or by populating different internal hyperfine 
states for cold atoms. Given the recent breakthroughs in 
the control and measurement of single optically trapped 
atoms [IHl [2DH^ . the experimental verification of the 
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coherent collective behavior of many particles, which we 
will show to heavily contrast the familiar established sta- 
tistical effects, is in reach. 



II. FRAMEWORK 

Single-particle evolution 

For our focus on many-particle interference, we exclude 
any interparticle interaction. The many-particle behav- 
ior can then be inferred from the single-particle time- 
evolution by appropriately summing the emerging many- 
particle amplitudes or probabilities. Under the action of 
the single-particle time-evolution operator U, an input 
state 10") evolves into a superposition of output states 

n 

fc=i 

where a (b) refers to input (output) modes, and j (k) is 
the respective mode number; the first (second) index of 
any scattering matrix refers to the input (output) mode. 
The probability 

P3,k = |C/„fcP, (2) 

for a particle in the jth input mode to reach the kih. 
output mode contains all occurring double-slit like single- 
particle interference. For identical particles, it is useful 
to work in second quantization: The creation operators 
a}j and for particles in the input- and output-modes, 
respectively, inherit the relationship ([T]), 

n 

«1->EMI (3) 

k=l 

which is valid for, both, fermions and bosons. The 
behavior of many-particle states is governed by the 
(anti)commutation relations for the creation and anni- 
hilation operators for (fermions) bosons. 

Particle arrangements 

Initially, N particles are prepared in the n input 
modes, with rj particles in the j-th mode, i.e. X]j=i ''j ~ 
N . We denote such arrangement of particles in the modes 
by a unique input mode occupation list r = (ri, . . . , r„) of 
length n. For distinguishable particles, the specification 
of the particle arrangement leaves the freedom to dis- 
tribute the (labeled) particles among the input modes in 
several different ways. Without restriction of generality, 
we choose the initial state 

i*'")d = ®u «4i \r,)) , (4) 



in the first-quantization formalism. For bosons and 
fermions, the specification of a mode occupation list r 
specifies the initial quantum state, 

i*")f/b = n M=r ' 

in the second-quantization formalism, where |0) denotes 
the vacuum. 

After time-evolution according to ^ and ([3|, respec- 
tively, the number of particles in each output mode is 
measured. Events are characterized by the corresponding 
particle arrangement defined by the output mode occupa- 
tion list s = (si, S2, ■ • ■ , Sn), where again — ^■ 

It is convenient to define for each (input or output) 
arrangement q an alternative notation d{q), the mode as- 
signment list. The list is of length N, with entries that 
specify each particle's provenience or destination. It is 
constructed by concatenating qj times the mode number 
j- 

d{q) - (i) = (1^_^, 2__J, . . . , rv^)(6) 

91 92 9t. 

The relationship between mode occupation- and mode 
assignment lists is exemplified in Fig. [T] 

Particle arrangements with cyclic symmetry will play 
an important role in our subsequent treatment. For any 
integer m that divides n, we define an arrangement q of 
N particles to be m-periodic if it consists ofp — n/m rep- 
etitions of a pattern k of length m. The mode occupation 
list thus reads 

q = {k\,k2, ■ . ■ , k„i, fci , . . . , k„i, . . . , fci , . . . , fcm), (7) 

" V ' 

p—n/rn 

while the mode assignment list d{q) satisfies 

Vj : dj+N/p{q) =dj{q)+m , (8) 

where we identify djv+j = dj and dj = dj + n. For ex- 
ample, the arrangement q — (2,1,0,5,0,2) exhibits no 
periodicity, and m = n = 6,p=l; the strongest symme- 
try is exhibited by an arrangement with Vj : qj = ri, as, 
e.g., q = (4, 4, 4, 4, 4), for which p = n,m — 1. 

Transition probabilities 

Signatures of many-particle interference appear in the 
transition probability, 

^D/B/F(r,s;C/), (9) 

for an input arrangement f, to an output arrangement 
s, given the single-particle evolution U . Bosons (B) and 
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FIG. 1: Many-particle paths and single-particle analogy, (a) The mode occupation hst of the initial (upper) arrangement 
reads r = (2,0,0, 1), the equivalent mode assignment list is d{r) = (1, 1,4), see Eq. ([g]). The final arrangement is characterized 
by s* = (1,1,1,0) and d{s) = (1,2,3). Three physically distinct many-particle paths, each corresponding to a permutation 
a of d{s), connect the arrangements. The many-particle (path) amplitudes are the corresponding products of single-particle 
amplitudes, given here below the respective paths. They are summed for all possible permutations in Eqs. (10 1 and (141 to 



give the full transition probability and amplitude, respectively, (b) Single-particle analogy: A single particle can pass through 
three slits, which corresponds to three distinct complex amplitudes. Again, the initial and final state are connected through 
three single-particle paths. 



fermions (F) can thus be compared to distinguishable 
particles (D). 

For the latter, the probability to find the final arrange- 
ment s is obtained combinatorially by taking into account 
all possibilities to distribute the particles among the out- 
put modes, given the single-particle probabilities pj_k in 
Eq. (|2|: 



TV 



Puir,s;U) 



E n 



(10) 



where sgnB((7) = 1 and sgnp{a) = sgn(cr) allows for the 
fermionic anti-commutation-relation. The transition am- 
plitudes for bosons (fermions) can be re-written as the 
permanent (determinant) of M: 



PBir,s;U) 



1 



|perm(Af)|' 



PFif,s;U) = |det(M)|2, 



in close analogy to (12 1. 



(14) 
(15) 



where S^-^^ denotes the permutations of the output mode 
assignment list d{s). We define the N x N matrix 



such that 



PD(r,s;C/) 



fperm(|Af|2), 



(11) 



(12) 



where the absolute-square is understood to be taken 
component-wise and perm(|Mp) den otes the permanent 
of |Mp 1231 mi- Each summand in (10) represents one 



way to distribute the particles among the output modes 
with Sfc particles in the fcth mode. Each possibility corre- 
sponds to a many-particle path as illustrated in Fig.jlja). 

For indistinguishable particles, all many-particle paths 
contribute coherently to the final state, and their ampli- 
tudes need to be summed, such that 



PB/F(r,s;[/) 



(13) 



N 



E sgnB/F(o-)n^rf.«,'^0-) 



Single-particle analogy 

An analogy between coherent and incoherent many- 
particle processes with the corresponding single-particle 
phenomena is suggestive: Consider a single particle that 
can pass through A^! distinct slits to fall onto a chosen 
point on a screen, as sketched in Fig. [l] If the time 
evolution is coherent, i.e. if the path taken by the parti- 
cle is not observed, the amplitudes of the A^! paths are 



summed, just like in Eq. ( 13 ). The probability to observe 
the particle reflects changes of the relative phases of the 
superimposed path amplitudes. Since the phases that 
are accumulated between each slit and the observation 
point can be, in principle, adjusted independently, fully 
destructive or strongly constructive interference can be 
induced. If the time-evolution, instead, occurs in an inco- 
herent way, because, e.g.^ the path information is leaked 
to the environment, the respective path probabilities need 
to be added, and the dependence on the accumulated 
phases vanishes, similar to Eq. (10 1. 



In the many-particle situation, we deal with a total of 
°'- ' n-1 J {n-l)\Nl ^ ' 



arr arr 



distinct particle arrangements for bosons or distinguish- 
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able particles, and 



NU 



Ny. 



(17) 



fermionic arrangements, i.e. the number of possible 
counting events rapidly grows with the number of par- 
ticles and modes. Again, the variation of the phases that 
are accumulated by the particles can be monitored by 
the many-particle event probabilities Pb/f given in (13 1. 
In contrast to the single-particle scenario, the up toTV! 
different many-particle amplitudes that enter in ( 13 1 are 



not independently adjustable, since they are all given by 
products of single-particle matrix elements Uj^k- There- 
fore, the behavior of many-particle scattering systems of 
moderate size {N k, uk, 10) already presents a challenge, 
which is aggravated by the large number of events ( 16 ) 
and (IT7|. 



The aforementioned scaling argument also immedi- 
ately exposes the main challenge for the exploitation of 
many-particle interference: For fermions, the Pauli prin- 



ciple implies N < n, and Af , given in ( 11 1, is a submatrix 



of U. Transitions are thus suppressed when M is singular, 
according to Eq. ( 15 ) that relates the transition ampli- 



tude to the determinant of M. For bosons, the perma- 
nent of M governs the behavior, for which no analogous 
criterion for its vanishing exists |25) . For distinguish- 
able particles, we face the permanent of a positive matrix 



in Eq. (12), which can be approximated efficiently [26] . 



The design of a setup that exhibits strong and controlled 
many-particle interference is thus in general much more 
involved than in the single-particle case. 



III. SUPPRESSION LAW 

A systematic assessment of many-particle interference 
becomes possible by imposing symmetries on the scatter- 
ing setup. We therefore focus on the Bell multiport beam 
splitter [27] , which is described by an unbiased scattering 
matrix in the sense that all single-particle probabilities 
are equal, |J7fc,zp = Pk,i = ^/n. The scattering matrix 
describes a discrete Fourier transformation, 
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(18) 



The phase that a particle acquires thus depends on the 
input and on the output mode, and can only assume a 
multiple of ir/n. In this setup, we expect a large vis- 
ibility in the sense that fermionic/bosonic probabilities 
strongly differ from their counterparts with distinguish- 
able particles, since many amplitudes of equal modulus, 
but of different phase, are added in Eq. (13). For distin- 



guishable particles, Eq. ( 10 ) gives a purely combinatorial 
expression: 



Pu{r,s;U'n = 



m 



.1' 



(19) 



'■3=1 . 

i.e. events occur with probabilities according to a multi- 
nomial distribution [23', which generalizes the binomial 
distribution found for a two-mode beam-splitter [28 . 

For bosons and fermions, the evaluation of the transi- 



tion amplitude ( 13 ) is, in general, not significantly sim- 

than for the 



pier for the highly symmetric matrix ( |18[ ) 

general case [29]. As a result of (U^°^) 
can, however, exploit an input-output symmetry: 



we 



^b/f( 



r,s;C/r') = ^B/F(s,r-;[/,f") 



(20) 



Additionally, arrangements that can be related to each 
other by a cyclic permutation or by reversing the mode 
order are equivalent, i.e. the application of these trans- 
formations on the initial or final state does not change 
any event probability. 

Combining the symmetry properties of the Fourier ma- 
trix (18) for m-periodic {i.e. cyclicly symmetric) initial 
or final arrangements, as introduced in Eqs. ([7|, ([8|, we 
can formulate a sufficient criterion for the occurrence of 
fully destructive interference and thus for the strict sup- 
pression of a transition r -f-?- s , which very considerably 
generalizes a results for N bosons that arc prepared in 
Tc = (ljlj---jl) and transmitted through an iV-port 
beam splitter [TBI ISD] ■ 



For bosons, given an m-periodic initial state r, final 
states s are suppressed when the sum of the elements of 
their mode assignment list, d{s), multiplied by the period 
length of the initial state, m, is not dividable by n, i.e. 



mod 



N 



dj{s),n ^ =^ PBif,s;U'n=0. 



(21) 



For fermions, the anti-commutation relation leads to a behavior that depends on the parity of the number of 
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particles: 



N 



N odd, or N/p even: mod m^<ij (s),n 7^ => Ppl^, s; C/^°") = 0, 



N 



N even and N/p odd: mod | m^dj(s),n ] ^ - => PF(f, s; ?7^°") = 0, 



r 



(22a) 
(22b) 



where n is even when N is even and N/p is odd. The 
proof for the suppression law is given in the Methods sec- 
tion. Intuitively speaking, Eqs. (pTj) and (22) formalize 



the observation that amplitudes of equal modulus but 
different phase annihilate each other when they are dis- 
tributed in an equally spaced way on a circle around the 
origin in the complex plane. 



The suppression laws (21), (22) circumvent the inher- 
ent complexity of expression (13): Even for very large 



particle numbers, the suppression of an event can be pre- 
dicted easily, whereas the computation of the permanent 
( 14 ) is computationally hard and thus practically unfea- 
sible for combinations of large N and n. For fermions, the 
determinant (15) is computationally less expensive than 



the permanent, but the suppression law (22) still offers 



a sig nificant speedup: The evaluation of the determi nant 
) by the LU decomposition scales with A'"'^, while (22) 
es linearly in 
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sea 



N. The suppression law also encom- 
passes several other criteria for the suppression of bosonic 
events in the literature as special cases [121 HH ISHHS^ . It 
does not give, however, a necessary criterion: It is possi- 
ble to arrange amplitudes in the complex plane such that 
they cancel each other, while they do not lie on a circle. 

Formally, the suppression of events given a cyclically 
symmetric initial state can be interpreted as the mani- 
festation of periodicity of the initial state in its Fourier 
transformation, given by Eq. (18). Very distinct initial 



states that share the same periodicity exhibit the same 
strict suppression among the final states. This reflects 
their shared periodicity, whereas the information on the 
exact constitution of the arrangement is contained in the 
unsuppressed states. The single-particle transformation 
that corresponds to the Fourier matrix can also be seen 
as the transformation to the quasi-momentum basis |33j , 
and the quasi-momentum-distribution can be probed |34| 
to characterize the phase of a gas. 



five many-particle interference thus leads to a quan- 
tum enhancement larger than unity (reddish colors), 
whereas a value smaller than unity (blueish colors) in- 
dicates destructive interference. Fully suppressed events 
are marked in black (green) when they are (not) predicted 
by the suppression law (21). 



The arrangements are ordered according to their occu- 
pancy, such that final (initial) arrangements with many 
particles in few modes are found on the left (upper) part 
of the plot. Arrangements which can be related to each 
other via cyclic permutation or inverse are omitted (they 
lead to identical transition probabilities), such that 50 ar- 
rangements remain. Initial states with a period (marked 
with a red arrow) lead to many fully suppressed tran- 
sitions due to the suppression law, whereas the remain- 
ing unsuppressed events are, consequently, typically en- 
hanced. 

A coarse-grained trend emerges: Events with few, 
highly occupied modes (on the left side) tend to be en- 
hanced, events for which the particles are well distributed 
over the modes (on the right side) are rather suppressed. 
A closer look reveals that this average trend has many ex- 
ceptions: Suppressed transitions can be found within the 
events with many particles in few modes, and enhanced 
transitions with many occupied modes appear; fully de- 
structive interference occurs within the highly and within 
the sparsely occupied arrangements. In general, the be- 
havior of an event is sensitive to the exact state prepa- 
ration and on the final state configuration: Shuffling one 
particle into a neighboring mode often turns constructive 
interference into destructive, and vice-versa. Similarly, 
the interference pattern shown in Fig. [2] strongly depends 
on the phases of the matrix ( 18 ). Phase variations imme- 



diately impact on all transition amplitudes besides singu- 
lar cases in which only one many-particle path is possible, 
such as for the initial or final state (6, 0, 0, 0, 0, 0). 



Scattering of bosons 



Mixed initial states and the emergence of bosonic behavior 



With a scenario of = 6 bosons in an n = 
6 mode setup, we exemplify the impact of many- 
particle interference and of the suppression law ( [2T| ). 
An overview of the system is given in Fig. [2] where 
the quantum enhancement, i.e. the bosonic probability 



PBir, s; Uq°^) (Eq. ( |Al| ) divided by the classical proba- 



bility PD(r, s;C/|'°" )(Eq. (19)), is displayed. Construc- 



The interference effects in Fig. [2] are sensitive to the 
loss of indistinguishability of the particles. When the lat- 
ter is not ensured anymore, their interference capability 
is jeopardized, such that the resulting transition proba- 
bilities eventually approach the ones for distinguishable 
particles [55] . 

The interference pattern also depends on the coherence 
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FIG. 2: Many-boson interference in a setup with A'^ = 6 bosons and 



6 modes. The color coding ind icates 



tiie quotient of the probability for bosons, PB{r,s;U °"), Eq. (All, to the classical probability PD{r,s\U ^ ), Eq. (191, as a 
function of the input state r and of the output state s. The input- and output configurations are arranged in the same order. 
Black fields denote transitions that are suppressed due to the suppression law (|21l), while green fields represent suppressed 



events which are not predicted by the law. One can identify the periodic initial and final states (marked with red arrows) as 
the black horizontal and vertical structures. For the initial states r = (0,0,3,0,0,3) and r = (0,1,2,0,1,2), the same final 
arrangements s are suppressed, since the suppression law only depends on the period length of the initial state. As an exception, 
the transition r = (0, 1, 2, 0, 1, 2) to s*= (0, 2, 0, 2, 0, 2) is suppressed, whereas for f = (0, 0, 3, 0, 0, 3), it is not. This is rooted in 
the period of s*= (0,2,0,2,0,2), m = 2, and 2- dj ((0, 1, 2, 0, 1, 2)) = 50 7^ mod 6, i.e. the transition is suppressed owing 
to the reverse of the law, which is obtained by exchanging the input and output states. The behavior of r = (0, 0, 3, 0, 0, 3) 
and r = (0,1,2,0,1,2) for the remaining, unsuppressed transitions differs strongly. The first line (marked by a blue arrow) 
corresponds to the initial state r = (0, 0, 0, 0, 0, 6), which does not exhibit any interference effects, since all particles are initially 
in the same single-particle state and no competing many-particle paths with different phase arise. A hierarchy of initial states 
can be observed for the final state s = (0, 0, 0, 0, 0, 6) (outmost left column): The more spread out the particles are in the input 
modes, the stronger is the enhancement of the final state due to bosonic bunching. 



of the many-particle wavefunction, i.e. on the purity of 
the initial state and on the coherence of the time evolu- 
tion. In order to see this, we consider a state in which the 
individual particles are still fully indistinguishable, while 
they are prepared in a mixture of all possible nonequiv- 
alent arrangements, {i.e. arrangements that cannot be 



related to each other via cyclic symmetry), 
1 



Pn 



arr 



(23) 



where N^^^ is the total number of arrangements given by 
Eq. (16). The probability for a final arrangement s for 
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Pmix amounts to 

f (Pmix,s) 



(24) 



i.e. to the transition probability averaged over the 
nonequivalent initial states r. The interference effects 
that govern the behavior for each pure initial state r are 
averaged out for pmix, and all final arrangements become 
approximately equally probable; we expect 



1 



7VB • 

arr 



(25) 



We verify this statistical argument by compar- 
ing the event probability for distinguishable particles, 
Pd(s, ^^^°"), the probability for bosons prepared in the 
arrangement fc = (1,1,1,1,1,1), the estimated prob- 
ability and the probability given the mixed initial 
state yOmix in Fig. [3] When the particles are prepared in 
fc, i.e. in a pure state, the system exhibits interference, 
and the resulting event probabilities strongly differ from 
the combinatorial value ( 19 ) for distinguishable particles. 



In comparison to the latter, events with many particles 
per mode (left side) are enhanced, events with many oc- 
cupied modes (right side) are rather suppressed in the 
equiprobable distribution Pe given in Eq. (25 1. Finally, 
the probability distribution for the state Pmix, Eq- (24 1, 



matches Pe very well, which confirms the approximation 
(25 1. Due to the mixed preparation of the particles, no 
interference pattern can crystallize for Pmix, and only an 
average, bosonic behavior remains. 

In other words, we recover the well-known statistical 
behavior of bosons when the purity of the initial state is 
destroyed, while the behavior of distinguishable particles 
emerges when the indistinguishability of the particles is 
lost. 




20 30 
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FIG. 3: Behavior of pure and mixed states of A'^ = 6 
bosons that are scattered in an n — 6 mode setup. 

We show the probability for final arrangements s, arranged 
in the same order as in Fig. [2] for distinguishable particles, 
Pd (blue circles, dotted line, Eq. (19l), for the equiprob- 
able distribution Pe = l/A'^^^ (black squares, solid line, 
Eq. (25|), for the mixed state of bosons pmix (Eq. (23l) 
to which all possible initial arrangements contribute with 



the same weight, i.e. {PB{f,s;U„°^))f (red diamonds, dash- 



dotted line, Eq. (24l), and for bosons that are prepared in the 
initial state fc — (1, 1, 1, 1, 1, 1), Pb(t^c, s; (7|^°") (green trian- 
gles, dotted line). Since we group each arrangement s* together 
with its cyclically permuted and mirrored counterparts, Pe (s) 
is not a constant, but it reflects the multiplicity of equivalent 
arrangements. The trend to favor arrangements with large 
populations is visible by comparing the probability for dis- 
tinguishable particles Pd to the estimate Pe, which favors 
multiply occupied states. The exact calculation for bosons in 
fc, however, does not exhibit this trend at all, and the pic- 
ture is dominated by many-particle interferences. Only when 
the mean over the initial states is performed, as in Eq. ( 23 1 , 



we lose many-particle interference and recover a clear bosonic 
bunching tendency, reflected by the values of pmix. 



Boson-fermion comparison 

When we restrict ourselves to initial and final states 
with at most one particle per mode, which we name Pauli 
states in the following, we can compare distinguishable 
particles, bosons, and fermions. 

In any setup with two particles, exactly two two- 
particle paths contribute to the total transition ampli- 
tude between any initial and any final Pauli state, which 
is reflected by the two terms in the sum (13 1. The two 



amplitudes are summed for bosons and subtracted for 
fermions, which leads to an antipodal behavior. For three 
or more particles, however, more paths contribute to the 
amplitude in (13), and no dichotomy is observed. 

Indeed, when we focus on the Fourier matrix (18), 



the very same transitions are predicted to be suppressed 
for an odd number of bosons and fermions, according 
to Eq. (21) and (22a). This emphasizes that our cri- 



(anti)commutativity of (fermionic) bosonic operators, 
but rather on the coherent superposition of many-particle 
paths. 

At first sight, it seems that the suppression law for even 
particle numbers and odd N/p (22b) predicts an anti- 
podal behavior of bosons and fermions. For exactly two 
particles, any arrangement with cyclic symmetry has in- 
deed m = n/2,p — 2, and transitions that are suppressed 
for bosons are enhanced for fermions, and vice-versa. In 
general, transitions with 



Q := mod m^dj(s),n g{0,7i/2}. 



(26) 



terion for destructive interference does not rely on the 



are only suppressed for one species, and not necessar- 
ily for the other. For large n, however, most transitions 
naturally lead to values of Q different from and n/2. 
Consequently, fermions and bosons also share many sup- 
pressed events for even particle numbers, unless N = 2. 
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FIG. 4: Many-boson (a) and many-fermion (b) interference of A'^ = 4 particles in a n = 12-mode setup. The 

color code reflects the quotient of the event probability for bosons (a) and fermions (b), divided by the average probability to 
find a Pauli state. The arrangements are given by their mode assignment list, the order of the final arrangements is the same 
as for the initial arrangements. Three out of the 29 inequivalent states possess a cyclic symmetry. The arrangements with 
elementary mode occupation lists k = (0,0,0,0, 1, 1) and k — (0,0,0, 1,0, 1) have period length m = 6, such that the same 
transitions are predicted to be suppressed for bosons and fermi ons; Eqs. (21 1 and (22a I apply. In contrast, the arrangement 



based on k — (0, 0, 1) has period length m — 3. Since N/p = 1, (22b I applies 



Mean bosonic and fermiontc behavior 

Despite the similarities between fermionic and bosonic 
interference, a mean trend for the probabihty to yield a 
Pauli state appears: Since events other than Pauli states 
are strictly forbidden for fermions, the average proba- 
bility for Pauli states is enhanced with respect to dis- 
tinguishable particles. For bosons, events with higher 
occupation are favored, and Pauli states are rather sup- 
pressed. 

In a Bell multiport with n = 12 modes and iV = 4 par- 
ticles, we find this tendency: The probability for a Pauli 
state amounts to 1/864 for distinguishable particles, ac- 
cording to Eq. ( |19[ ). The average fermionic probability 
is enhanced, it amou nts to the inverse of the number of 
available states, Eq. (jl7|, (iVf„)~^ = 1/495. For bosons, 
events with multiple occupation are privileged, the aver- 
age probability for a Pauli state becomes w 7.50 • 10"'', 
which is close to the inverse of the number of accessible 
states, {N^,,y^ = 1/1365 « 7.33 • IQ-*. 

From this general trend, however, no systematic state- 
ment can be inferred for the individual transitions, which 
are displayed in Fig. [4} No anti-correlation between 



fermionic and bosonic behavior is visible, and the correla- 
tion coefficient between the enhancement or suppression 
of fcrmionic/bosonic events amounts to only -0.05. The 
enhancement of a fermionic transition does not imply the 
suppression of the corresponding bosonic transition. 



Decreasing the particle density, i.e. increasing the 
number of modes n for a constant number of particles 
N, has no impact on the strength of interference ef- 
fects: Fully suppressed and highly enhanced transitions 
remain, and so does the general impact of many-particle 
interference. In turn, the difference between the average 
bosonic/ferminionic behavior fades away when the num- 
ber of modes is increased: The quotient of bosoni c to 
fermionic accessible states, N^^.^. / N^^.^. , given by Eqs. (16) 



and (17 1, approaches unity for n — > oo, and the average 



statistics then resembles the combinatorial distribution 
of distinguishable particles. This reflects the absence of 
the Fermi pressure [2] and of bosonic behavior for low- 
density systems. 
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FIG. 5: Many-boson and many-fermion scattering 
of pure and mixed states. The probabilities for fi- 
nal arrangements s are arranged in the same order as in 
Fig. |4] Bosons (green triangles, dotted line) are compared 
to Fermions (brown circles, dashed line) for the initial state 
d{fp) — (1,4,7,10) that corresponds to the very last line in 
Fig. |4] An average is taken over all nonequivalent initial 
states, which results in the black squares for fermions and 
red diamonds for bosons. It approaches the inverse of the to- 
tal number of states for fermions (blue solid line) and bosons 
(red solid line), respectively. 



Mixed initial states 

Just like for the case of six bosons discussed above, the 
transition amplitudes for fermions and bosons that are 
prepared in Pauli states depend on the initial pure prepa- 
ration of the particles: Given a fully mixed state analo- 
gous to (23), but restricted to PauH states, the resulting 



event probability approaches a constant value that does 
not depend strongly on the final arrangement. Since the 
indistinguishability of the particles is not jeopardized by 
the mixedness of the initial state, the Pauli principle pre- 
vails for fermions, and so does the average privilege of 
multiply occupied states for bosons. The probability for 
any final state then approaches the inverse of the num- 

ber of states, ( TVa/r ) , which were given in Eqs. ( 16 

[It]). The behavior of bosons and fermions is compared in 
Fig. [5j where the loss of interference due to the mixedness 
of the initial state is apparent. 

In analogy to the bosonic case, the loss of the purity of 
an initial state of fermions leads to the loss of the inter- 
ference pattern, while the quantum statistical behavior 
persists. Only when the indistinguishability of fermions 
is jeopardized, also the impact of the Pauli principle fades 
away. 



IV. CONCLUSIONS AND OUTLOOK 

The coherent many-particle evolution of non- 
interacting identical particles is governed by many- 



particle interference, with consequences that go far 
beyond bosonic (fermionic) (anti)bunching encountered 
in incoherent environments. The latter phenomena can 
be understood from the postulate that assigns every 
micro-state the same realization probability [2 when the 
specific constraints to bosons and fermions are respected. 
The (anti)symmetrization, however, establishes also a 
many-particle coherence property that leads to the 
encountered interference effects, which seldom reproduce 
the familiar statistical behavior. 

Although any many-particle interference setup can be 
treated by Eq. (13 1, the analysis becomes rather tedious 



when the number of particles and modes are not both 
small, which is due to the computational complexity in- 
herent in Eqs. (14) and (15), and also due to the large 
number of accessible states, Eqs. ( 16 ), ( 17 ). For Bell mul- 



tiport beam splitters, we can circumvent these difficulties 
by exploiting the available symmetries, which allows the 
systematic confrontation of many-boson to many-fermion 
scattering. Since for any number of particles iV and any 
number of modes n with a non-trivial greatest common 
divisor, i.e. GCD{N, n) > 1, non-trivial periodic states of 
the form ([7| can be found, numerous applications of the 
suppression laws are possible. We thus have a construc- 
tive recipe for many-particle large- visibility setups at our 
hands, since events that exhibit destructive interference 
can be found for arbitrary system sizes, whereas the ex- 



plicit evaluation of ( 13 1 becomes prohibitive already for 
moderate particle and mode numbers. This provides a 
powerful characterization toolbox for the indistinguisha- 
bility of many particles, for the purity of their initial 
state preparation and for the many-particle coherence of 
the time-evolution - which are all imperative properties 
for quantum technologies. 

The interference patterns shown in Figs. 2][5 are spe- 
cific to the Bell multiport, but also exemplary for any 
other scattering scenario. Our results suggest a forceful 
distinction between quantum statistical effects^ i.e. the 
(anti-)bunching of uncontrolled (fermions) bosons in a 
thermal state [8] or in the mixed state ( [23| , and many- 
particle interference^ i.e. the coherent conspiration of 
many-particle paths. The former are the consequence of 
the interplay of the kinematic constraints with the sta- 
tistical uncertainty, whereas the latter are jeopardized 
by such uncertainty. A many-particle path picture is 
also applicable in the incoherent case [51 [S]. However, 
one then only distills an average, incoherent statistical 
behavior that can also be understood from the equal a 
priori probability postulate [2] . 

To name a further difference, quantum statistical ef- 
fects fade away with decreasing particle density, in con- 
trast to many-particle interference. So far, the literature 
has mainly concentrated on the evolution of two-particle 
states [TU] and on two-point correlation functions [SJ [5] , 
for which these phenomena appear to be synonymous: 
The Hong-Ou-Mandel effect 10] is in ostensible agree- 
ment with the bunching behavior of thermal bosons [3], 
and it also refiects a strong boson- fermion dichotomy [llj . 
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As we have shown, this ostensibly intuitive picture needs 
to be abandoned immediately as soon as more than two 
particles are considered. 

The loss of purity in the initial state leads to a dete- 
rioration of interference effects and to a recovery of an 
average bosonic or fermionic behavior, while a loss of 
mutual particle indistinguishability leads to the combi- 
natorial behavior of distinguishable particles. In gen- 
eral, decoherence is thus a more diverse and complex 
phenomenon in the many-particle domain than in the 
single-particle realm. It remains to be studied how deco- 
herence processes affect many-particle interference: We 
speculate that different mechanisms deteriorate the in- 
terference pattern in different ways, such that one re- 
covers - in the limit of strong decoherence - either the 
behavior of fully distinguishable particles [35J or a 
bosonic/fermionic statistical behavior. 
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Appendix A: Proof of the suppression laws 



For bosons and fermions, by inserting the definition 
of the Fourier matrix (18), we can rewrite the transition 



amplitudes ( 14 1 and ( 15 ) as sums of amplitudes with unit 
modulus: 



where sgng((7) = 1 and sgnp(cr) = sgn(o'). 



sgnB/F(o-)exp ( i'^Yd„U){r) ■ dj{s) 



N 



(Al) 



With the definition 



Bosons 



N 



e(f, s, cr) = ^ d„(j)(r)dj(s), 



(A2) 



the probability P^{r,s;U^°^'^) in (All becomes propor- 
tional to 



^ exp U — e(f,s,cr) j 



(A3) 



In other words, the natural number 8(r, s, a) is the total 
phase that is acquired by the many-particle wavefunction 
(in multiples of 27r/n), when one specific many-particle 
path defined by cr, from the initial state f to the final 
state s, is realized. Since 0(r, s, ct) is a natural number, 
the sum in ( A3 ) contains only n-th roots of unity. 



PB(r,s;f/F°") cx 



ji-i 

E 

fe=0 



Cfce 



(A4) 



with Cfc e N, X]fc=o "^fc = ^'i ^-nd where the value of Cfe 
corresponds to the number of permutations a for which 
mod {Q{r,s,a),n) = k P]: 

Uk{r,s) — {a\Q{f, s, a) — k mod n} , 

Ck = \uk\, (A5) 



where x = y mod k means that there is a integer number 
I such that X = I ■ k + y. Hence, the sets Uk group all 
many-particle paths that acquire the same phase, 27rfc/n. 
Since all moduli of involved probabilities are equal, these 
paths possess the same amplitude. It is helpful to define 



/ N \ 

Q{m, s) = mod j m di{s),n ] 



and an operation F on the permutations cr. 



r(a)(fc)= (a(fc) + ^ 



mod 



(A6) 



(A7) 



which shifts the permutation cr by the number of par- 
ticles in each period repetition, N/p. Due to the m- 
periodicity of the initial state r, which implied ([8|, the 
value of the total acquired phase Q{f,s,a) acquires the 
constant Q{m, s) when the above transformation is ap- 
plied on a permutation a: 



N 



e(f,s,F(a)) = 5]d,(s)dr(.0))(rO. (A8) 
Using ^ and ( [AT] ), we find 

N 

e(f,s,F(CT)) = Ydj(s)d„(^j){f 



N 



9 (f, s, cr) -I- Q{m, s) mod n. (A9) 
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Thus, if Q{m, s) 0, the repeated application of F gives 
us a bijection between all pairs of Ub+a-Q{m,s) 

for b e 

{1, .., n} and a G {0, 1, .., ti — 1}. Therefore, we have 

V6 e {0,...71- l},Va e IN : Ci,+aQ{rn,s) = Cb, (AlO) 

Since only index values mod n are relevant, all equalities 



that can be inferred from (AlO) are contained by the 
following n equalities: 



V6 e <iO,. 
Va e <!0, . 



n ■ Q{m, s) 
LCM {Q{m,s),n) 
LCM (Q(m,s),n) 
Q{m, s) 



- 1 



1 



(All) 



where LCM(x, y) denotes the least common multiple of 
X and y. Since 



LCM(x,y)GCD(a;,?/) =x-y, 



(A12) 



a. Odd particle number N , or even N/p 

For odd particle numbers N, the application of F on 
a permutation a does not change its parity: Any cyclic 
permutation on a set of odd size (here the set of indices 
{1, . . . , N}, of size N) is itself an even permutation be- 
cause it is composed of the even number — 1 of (odd) 
elementary transpositions. 

The parity of a permutation remains also unchanged 
under the application of F when the number of particles 
N is even and N/p is also even. The application of F on a 
permutation corresponds to an N/p-io\d cyclic shift that 
is composed by N/p (odd) elementary shifts. Its parity 
is therefore even, since N/p is even. 

In these two cases, we find, as for bosons but now in- 
dependently for even and odd permutations: 



b can effectively take GCD{Q{m, s), n) —: g distinct val- 
ues, and each Cb is equal to n/g — 1 other coefhcients 



Cb+aQ- By setting the summation index k in (A4) to 
k = b + a ■ Q, the sum ( A4| can be rewritten as 



PB(r,s;C/P™)oc 




Q{rn,s)-a 



(A,13) 



where we exploited explicitly (AlO), such that the sum 
factorizes in two parts. The sum over a is a truncated 
geometric series. 



3=0 



1 - a;' 



(AM) 



with I — n/g and x — e*^'^'^™'^. Since x is a n-th root of 
unity and x 1 (since we assumed Q{m, s) ^ 0), ( A14) 
vanishes, and so does (A13). 



Adaptation for fermions 



V5 e {0, . 



1} 



„E P 
^b 1 C^+Q 



(A17) 



The symmetry property of (AlO I is hence inherited in 



dependently by the even and the odd permutations. The 
formulation and the consequence of the suppression law 
consequently remain unchanged: Formally, it states that 
the even and odd part of the sum ( A15) both vanish. 



b. Even particle numbers N and odd N/p 

For even TV and odd N/p, the application of F onto a 
permutation a changes its signature, and we find 



Vfee {0, ...n- 1} 



-b+Q 



-b ■ 



(A18) 



The signatures of the permutations are thus inter- 
changed, and our suppression law in the formulation (21 ) 



is not valid anymore. By a case-by-case analysis, we can 
explore the consequences of (A18|: 



In order to adapt the suppression law for fermions, we 
need to include the signature of the resp ective permuta- 
tion in the sum of the amplitudes (Al). In analogy to 
(lAil), we find 



Pf (r, s; U 



Fou\ 



cx 



oi^O(f,s,cr) 



(A,15) 



where we split the sum ( Al ) into even permutations S^, 
and odd ones S^j- We now need to consider the sets of 
even (E) and odd (O) permutations with 8(r, s, a) — k, 



E{0) 



(f,s) 



(A16) 



{a\Q{f, s,cr) = b mod n, sgn(CT) = +1(— 1)}, 



separately. Their respective cardinality is denoted by 
= [wf^*^''!. We also need to infer the action of F on 
the parity of the permutations. 



• The condition mod{Q, n) = is sufficient for the 
full suppression of the respective transition. In this 



case, V6 



and each amplitude cf'e' " 



has an amplitude — c)j e*^ equal in magnitude, but 
opposite in sign. A transition between r and s that 
is not necessarily suppressed for bosons becomes so 
for fermions. 



n/2 is not neces- 



• A transition with mod{Q, n) 
sarily suppressed, since = c^^i, such that two 

amplitudes and c^^n do not cancel - as they 
do for bosons where they lie on opposite sides of 
the origin in the complex plane - but instead they 
enhance each other. 



For all other values of mod{Q,n), i.e. values that 
are neither zero nor n/2, the respective transi- 
tions are, again, necessarily suppressed, since ( A13) 
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turns into 



PF(r,s;C/^°")(x 



\r=0 

n/g-1 

a=0 



, (A19) 



where the sum over a can, again, be represented as 
truncated geometric series and vanishes since n is 
even by assumption. 
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